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АВ5ТКАСТ 


An investigation into the problems of economically pro- 
cessing sewage on board Naval ships resulted in the develop- 
ment of two computer simulations employing Monte Carlo 
analysis to describe the generation of sewage. Simulation 
one was based on a non-homogeneous Poisson process. For 
simulation two, an empirical distribution describing the 
arrival behavior of sewage to the holding/processing unit 
over a 24 hour period was applied to known data on sewage 
generation. 

Results of the two simulations were compatible with one 
another. Aside from pointing out a most feasible combina- 
tion of holding tank capacity, processor rate and processing 
Осу, the simulations also indicated that a revision of 
the Navy's design parameter for the daily per capita sewage 
generation rate was in order. The simulations were designed 
so that use by others with their own data would be an easy 


matter. 
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1. INTRODUCTION 


A. BACKGROUND 

During the time between the conception of this thesis 
and its completion, significant changes in both policy and 
attitude have been made by the Navy in the field of pollu- 
tion abatement. These alterations have resulted from in- 
‘creasing public pressure for the-ecological restoration of 
U.S. inland waters and subsequent congressional legislation 
which outlined future restrictions to all seagoing vessels 
and watercraft. More specifically, in August 1971, the 
Chief of Naval Operations stated "that the water pollution 
abatement program for Naval ships has as its long term goal 
the installation of two basic systems aboard ship: 

1. A Marine Sanitation Device (MSD) which processes 
sewage to an acceptable effluent; 

2. A Collection, Holding and Transfer (CHT) System 
which is designed to transfer all liquid wastes (except 
oil) ashore or to a barge." 

Late in 1971, however, the Senate issued further re- 
strictions to vessels to the extent that, by the year 1981, 
there shall be no toxic effluent discharge, whatsoever, 
from any ship into U.S. inland waters. As a result, Naval 
design efforts refocused primarily on the development of 
a Collection, Holding and Transfer System. This shift in 
emphasis was also aided by the fact that prototypes of 


Marine Sanitation Devices had suffered from low system 





reliability and excessive costs. "The cpresentacgc си 
of the Navy is the development of a sewage holding tank 
system with the capacity to retain 12 hours of sewage gen- 
erated by a full complement crew and with the capability of 
transferring such sewage to an external processing facility 
ashore. Such a system is to be installed in all ships re- 


maining in the fleet after 1977. 


B, STATEMENT OF THE PROBLEM 

It is the contention of the author that neither strictly 
holding sewage for later transfer, nor immediately process- 
ing the sewage internally as it arrives to the processing 
unit, is the most feasible economic solution to pollution 
abatement. A policy of holding sewage and transferring it 
to some other processing station would greatly hamper the 
I y's capability to act as a mobile self-sustained unit; 
likewise, the policy of rapid processing aboard would re- 
quire expensive units of high output and rellability. It 
15 believed that viable trade-offs can be realized through 
a combination of the two, employing a smaller, cheaper pro- 
cessor which could operate at a lower processing rate, and 
a smaller holding tank that could absorb sewage flow during 
peak periods of the day. 

A second contention, that is perhaps more germaine to 
the current Navy policy of a twelve hour holding system with 
subseguent transfer, is the belief that the presently ac- 


cepted daily per captia rate of 30 gal/man/day of sewage 





generation is over estimated [Ref. 5]. The parameter was 
determined from data gathered by the Marine Engineering 
Laboratory on the daily sewage generation rates per man on 
board four ships. By reasoning that "any treatment unit 
used aboard a Naval vessel must be designed to operate with 
a full complement 24 hours per day," the MEL Report [Ref. 
18] recommended that the rates measured be doubled for de- 
sign purposes. This resulted in a figure of 30 gal/man/day. 
Since the Original figures were obtained on a per capita 
basis, and there was little variation, it is felt by the 
author that more accurate results could be obtained by 
applying the unmodified data to a model which also considers 
the traffic pattern of personnel using the facilities over 
a period of time. 

It was the purpose of this study, then, to investigate 
the problem of sewage generation and processing aboard a 
Naval ship, and to develop a model which would best de- 
Scribe such a process. This model would be useful in aiding 
decisions on the selection of the most feasible treatment 
stem. 

The approach taken, as presented in this thesis, was to 
first formulate the problem analytically and to compare it 
with other research performed in related cases. A computer 
simulation was then developed, motivated by a similar study 
done for the air lines regarding sewage handling on board 
large jumbo jets [Ref. 9]. The simulation, its methodoligies 


and how it can be. used are described in Chapter III. The 





simulation provides a sensitive model which illustrates the 
trade-offs obtainable through any treatment policy tested. 
It also reinforces the contention that a revision of sewage 
generation rates, combined with closer analysis, can lead 
to savings in holding tank size. 

It should be noted that a simulation is no stronger than 
the basic assumptions made regarding its underlying struc- 
ture, and that such assumptions OS examined criti- 
cally. With this in mind, the simulations were developed 
using the best information available; however, they were 
also made flexible enough so that, should a designer have 
better arrival data or other policies to test, such informa- 
tion could be easily read into the program to give useful, 


real world results. 





ТТ. PROBLEM FORMULATION 


A. THE BASIC MODEL 
The overall process of generating sewage over a specific 
length of time, storing it and/or processing it can be de- 
scribed analytically. Further, it can be related to a 
group of similar problems that are generally classified as 
storage problems, examples of which are the collection and 
storage of water in a reservoir subject to some release 
rule, and inventory problems in which demand is random. 
Any such process in which either the input or output is 
random comes under the category of a stochastic process. 
Such is the case with the sewage storage and processing 
problem, in which the input, X(t), (see Figure 1) into a 
storage tank is random in nature. Here X(t) represents the 
amount of sewage generated by the ship's company at discrete 
points in time, t, corresponding to, say, hours or minutes. 


Two important characteristics of this input are 1) its time 


Et) 
— 
Input = 
Z (t) 
Tank V Processor 
Figure 1. Schematic of System. 
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dependent arrival behavior to the tank, and 2) its overall 
distribution. It is assumed that X(t) ls independent and 
identically distributed for different values of t. For this 
discussion, X(t) is assumed to be solely bodily waste; how- 
ever the model could also be expanded to include "hotel" 

and galley waste. 

2 (t), represents the.total .amount of sewage in the stor- 
age tank also at time t. Aside from being a function of the 
input random variable, X(t), it is also a function of the 
rate, a, at which the sewage is processed, and the particu- 
lar "release rule" or processing policy, y, that is applied. 
The processing policy is essentially a switching function 
which turns the processor on or off based on either a time 
plan (e.g. on at = Е of the day), a tank level 
plan (e.g. turn on when the level of sewage reaches a cer- 
tain height), or a combination of the two. 

At certain times of the day, the processing rate may 
exceed the input, X(t), while at other times it may lag it. 
ТЕ the processing rate is much lower than the input rate, 
Z(t) will behave like a sum of independent random variables. 
For such a case, the tank would be expected to overflow. 

If the processing rate is much higher than the input, then 
the tank may never fill at all, and Z(t) will behave more 
deterministically (or equal zero). One of the questions to 


be answered, then, is what is the optimal processing rate? 


Pr 





In general, Z(t) may be described as follows: 


K if Z(t-1)+X(t) -at>K 
Z(t-1)+X(t)-at, if O«Z(t-1)-«X(t)-at«K 

Z(t) = (1) 
Е), if a=0 
0 if Z(t-1)+X(t)-at<o 


where K is the capacity of the tank. 
This formulation includes two special cases. The first 


is the event that overflow occurs in which case: 


Z(t) = 


Ep. if Z(t-1)4X(t)-at«K 
(2) 


К if Z(t-1)+X(t)-at>K 
Here the overflow is 7(+-1)+Х (+) -а=-К. The second case is 
the event that processing empties the tank completely, so 


that: 


Z(t-1)+X(t)-at, 1£ 2(t-1)+Xx(t)-at>0 
Z(t) = (3) 


0 if Z(t-1)+X(t)-at<0 


DR OTHER RELATED STUDIES 

Equation (1) through (3) are not unique to the sewage 
generation and collection problem. Rather, they parallel 
those proposed by Moran, Gani, and others [Refs. 15, 8, 10] 
under the broader discipline of the stochastic theory of 
reservoirs or storage, an approach which has attracted con- 
siderable theoretical interest, but has seen less application. 

Initially the theory, as proposed БУ POP required the 
inflow, X(t), to have a known continuous probability dis- 


tribution and to be mutually independent of the inflow for 
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any other value of t (all of which have the same distribu- 
tion). The level in the tank, Z(t), is subject to some re- 
lease rule, j|, and the sequence 1Z(t+1), t>0} is a Markov 
process, with state transition probabilities that can be 
computed from the inflow probability distribution. Consid- 
erable effort has been expended in finding solutions to the 
“Steady »"state distribution of Z(t) for both continuous and 
discrete cases. 

These initial studies were extended by a number of 
authors. One such work by F. M. Lloyd [Ref. 12] showed that 
the theory need not be limited to inflows that were mutually 
independent, but that it could also be applied to serially 
correlated inputs with the results forming a bivariate 
Markov Chain. Ali Khan and Gani [Ref. 1] carried this study 
one step further to consider the time dependent solution 
to the same model. 

An excellent example of the application of the above 
theory which is closely related to the sewage problem is 
the work by Avi-Itzhak and Ben-Tuvia [Ref. 3]. Here А 
authors studied the problem of "how big a reservoir to build 
and what should the pumping capacity be, such that the cost 
per cubic meter of water is a minimum." In order to obtain 
their solution, however, the authors first had to assume an 
approximate distribution for the daily input to the reser- 
voir as well as know the costs involved in operating the 


system. 


T3 





In all of the above cases the underlying requirement is 
that the distribution of X(t) is known or can be approxi- 


mated by a sequence of discrete distributions. 
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ТТТ. MODELING THROUGH SIMULATION 


А. JUSTIFICATION OF SIMULATION 

The storage theory approach described above leads to 

perhaps the most precise solution to the sewage storage 
problem. ТЕ has, however, certain drawbacks which preclude 
its immediate implementation. The first of these is the 
lack of real world information and data with which to 
derive the necessary matrix of transition probabilities. 
Secondly, the mathematics involved are difficult and cum- 
bersome. This is especially so when one is attempting to 
solve the problem for different processing policies ("re- 
lease rules") and processing rates, each of which would 
require its own transition matrix. Thirdly, the сте Л 
of the distribution of X(t), the input variable, whose sta- 
tionarity over time is questionable, would require parti- 
tioning the process in order to arrive at a solution. Such 
a non-homogeneous process would only further complicate 

the analysis. Lastly, the problem of a reversing boundary 
brought about by abrupt changes in processing rate, accord- 
ing to the processing policy used, would add even another 
dimension to an already difficult problem. 

Taking the above into consideration, it was decided 
that a computer simulation involving Monte Carlo analysis 
was at least a justifiable, if not preferable, alternative. 
Most important to the justification of a simulation was the 


fact that it could be made amenable to апу а Ее га ао Е аи 


T5 





processing policies and processing rates. It is also adapt- 
able to adjustments in the input distribution if and when 
additional data on the arrival patterns of sewage inputs 

are obtained. These characteristics have been incorporated 


into the simulations developed in this study. 


B. DESCRIPTION OF THE SIMULATIONS DEVELOPED 

‘Pwo computer ‘simulations were -developed using FORTRAN IV 
for runs on the IBM 360/67 computer. The simulations dif- 
fered only in the method of Monte Carlo analysis used to 
generate the sewage input distribution, X(t), which was then 
applied to the basic mathematical model described in Chapter 
II (equations (1) - (3)). For simulation one, X(t) was de- 
termined by assuming a non-homogeneous Poisson process. For 
simulation two, an empirical distribution describing the 
arrival behavior of sewage to the tank over a 24 hour period 
was applied to known data on sewage generation. The tech- 
niques and mathematics involved are described in supsections 
2 and 3 following an overall description of the basic simu- 
lation (subsection 1). 

1. Basic Design Common to Both Simulations 
Aside from the methods used to generate observations 

of the sewage input, the two simulations were identical in 
the way the sewage was processed from the holding tank. 
(See general flow chart of the simulations, Figure 2.) 
First it was assumed that the holding tank used in the 


models had sufficient capacity so that overflow was not a 
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GENERATE A 
RANDOM VARIATE 
FOR X(T) 













INCREASE LEVEL, Z(T) 
IN TANK ACCORDING TO: 
Z (T)2Z (T-1)4X(T) -aT 









INCREMENT 
TIME CLOCK 
T=T+1 







ING POLICY. EXAMINE 

(Z(T),T) TO DETERMINE 

IF PROCESSING STOPS 
(a=0) 


N END OF WEEK? 


X 
СОМРОТЕ STATISTICS 


PRINT 
RESULIS 


CALL THE DESIRED PROCESS- 


Figure 2. General Flow Chart of the Simulations. 
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problem. This was done in order to get some idea of what 
the maximum sewage level (ZMAX) in the tank would be, ac- 
cording to the processing rate and policy used. This 
further simplified the basic model, so that the only equa- 


tions needed for the simulation were: 


MN -at, a>0 
Z(t) = (4) 
0 if Z2(t-1)+x(t) -at<0 

Both simulations involved discrete time steps; for 
model one, the unit of time was 5 minute increments; for 
model two, hourly increments. In both cases the total period 
of time simulated was one week, based on the assumption that 
it was desirable to have any sewage build-up during the week 
processed, and the tank emptied by the end of the weekend. 
A new cycle would then correspond to the start of a new work 
week. It should be noted that each simulation day started 
at 0500 real time. At each discrete time increment, a ran- 
dom observation of sewage input, X(t), was generated and 
the level, Z(t), was computed using equation (4) for a 
specific processing rate, a. 

Anyone of five processing policies, in the form of 
separate subroutines could be called independently and ap- 
plied to the model. The following five policies were used: 


1) Plan 1. Continuous holding with no processing 
whatsoever (a=0). 


2) Plan 2. Continuous processing at some rate, a, 
to be read in by “ehewpreguan iam: 
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3) Plan 3. Process continuously from ОБОО nea 
1600 each day; after 1600 stop processing when 
the tank level reaches 50 gal; continue pro- 
cessing at 0500 the next day. For weekends 
stop processing anytime the tank level reaches 
50 gal; continue processing when the tank 
level builds up to 500 gal. (This is a com- 
bination of a level policy, Plan 4, and a time 
policy, Plan 5, as described below. 


A) Plan 4. (A level control policy) Stop pro- 
cessing whenever the tank level goes below 
„some, prescribed .limit (50 gallons was used in 
the simulations - the value is read in by the 
programmer). Continue processing when the 
level exceeds an upper limit (1500 gallons 
was used). 

5) Plan 5. (А time control policy) Process con- 

tinuously for a fixed portion of each day. 
(For the simulations processing stopped at 
hour 16, 2100 real time and continued at the 
beginning of the next simulated day, 0500 real 
time.) 

Policies one and two illustrate the two extremes 
possible. The policy of continuous holding (Plan 1) was used 
to determine the amount of sewage that would build up over a 
certain period of time. (The special case of a twelve hour 
holding tank was tested using this policy.) The policy of 
continuous processing (Plan 2) was used to get an indication 
of the minimum tank levels that could be expected over any 
range of processing rates desired. 

The remaining three policies were developed as ex- 
amples of probable real world policies. They also served to 
reflect the response of the simulations to changes in policy, 
for which the simulations produced acceptable results. Of 


the three, policy five would be expected to be the least ex- 


pensive since it simply involves turning the processor off 
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once a day. Policy four, on the other hand, would require a 
level sensing device and some sort of automatic start/stop 
mechanism to the processor and would, therefore, entail a 
considerably larger cost. Policy three would only require 
manual operation and its relative cost would fall somewhere 
in between that of policies four and five. The point to be 
ee three is the fact that the smaller the 
time required for the processor to operate, the lower the 
overall running cost and the time available for system 
maintenance. 

It is emphasized that the simulations are in no way 
limited to just these policies and that any other policy may 
be tested simply by the addition of a subroutine similar to 
those used and a C by which to call it. Further 
details on the mechanics of using the computer simulations 
are contained in Appendix D. 

2. Input Distribution by Simulation One 

One of the most important problems to be solved was 
how to generate reasonably accurate observation of the 
arrivals of sewage to the holding tank. The two parameters 
that had to be considered were 1) the arrival rate and 2) 
the quantity of sewage per arrival. Although a study was 
conducted in 1964 [Ref. 19] to determine the quantity of 
Sewage generated per man per day on a Navallyessel то 
definitive statistics were obtained to describe the arrival 
process over the course of a 24 hour period. Owing to the 


rature of shipboard routine, which can be considered 


20 





reasonably standard throughout the Navy, it was evident 
that the arrival distribution would be multimodal over the 
period of the day. With such changes in arrival rates it 
was also clear that the process could not be considered 
Stationary for a 24 hour period and that some form of non- 
"homogeneous process would be required. 

.In a.related study by Gephart and Balachandran 
[Ref. 9] into the reliability mI of an airborne 
sewage system for large jet liners, data was collected on 
the arrival behavior of passengers to the restroom facili- 
ties. It was concluded by the authors that the arrival 
process could be considered Poisson over "non-overlapping 
and consecutive time segments." It appeared reasonable, 
then, that a non-homogeneous Poisson process could also be 
applied to the sewage generation problem on board a Naval 
ship. 

In its most general form, the process could also 


be considered a compound Poisson process represented by 


М (+) 
X(t) = X Yi (5) 
i=1 
where {N(t), t20} is Poisson process and {Y,, n=1,2...) is 


a family of independent, identically distributed random 

variables representing the quantity of sewage per arrival. 
If a distribution for the body waste generated per capita 
per arrival were known, this could easily be incorporated 


into the simulation. For the purpose of simplicity, 
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however, the quantity of sewage per arrival was taken to be 
a constant and the compound property was disregarded. Based 
on studies [Ref. 19] which determined the quantity of water 
per flush to be 4.5 gallons, a rather conservative с 
of 4.9 gallons per arrival was used. 

Оро Тл! problem, then, was to determine the best 
«estimate. Of the..arrival .behavior at head facilities over a 
24 hour period. It was decided that the day could logically 
be partitioned into three periods: the first being the 
period between 0500 and 1600 when the majority of the crew 
is on board for the working day; the second being the period 
between "liberty call" and "taps" (1600-2200) and the last 
including the time from 2200 until 0500. The intensity 
function A(t) which describes the arrival rates could then 


be written for weekday periods as: 


À) Е. 
д (6) = Nat Deci (6) 
А», Eust 


and assuming that the majority of the crew is on liberty 
over the weekend, then a reasonable intensity function for 
weekends might be 


а (7) 


спете х) ала A, are identical fon БОЕВ ЕЕ вое ао 
Thus the mean valve: wre tone ce 
t 


m(t) ={ A(s)ds = AititA, (t,-t))+A, (t -t2) (8) 
О 


p 





t 
men) = а т = А (9) 
О 
and it can be shown that 


| e-m(t) [m(t)]" 
P{N(t)=n} = , n»0 (10) 
n! 





This then describes a non-homogeneous Poisson pro- 
cess. The values for each Л; were determined by comparing 
the mean value function m(t), multiplied by 4.9 gallons of 
sewage/arrival, to the expected daily and weekly sewage 
build-up. The details of this are contained in Appendix A. 


The final values used for each Àj are presented in Figure 3. 


Arrival/5 min 





0500 1600 2200 0500 ps 


Figure 3. Arrival Rates for Simulation One. 

The simulation involved generating Poisson variates 
for each 5 minute increments of simulated time, using the 
above parameters during the perlods Indica ceo TICA e 
of arrivals was then converted to a quantity of sewage X(t) 


according to; 
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X tjus NAE) 4.9 gal/arrival CI) 
giving the necessary random input for X(t) which is then 
applied to the basic model as described in section 1 above. 

3. Input Distribution as Described by Simulation Two 

Although the reasoning behind Simulation One was 
considered sound, it was decided that a second simulation 
could be developed through a different approach. Not only 
would it take full advantage of all known data, but it would 
provide an experiment by which the two simulations could be 
weighed against one another. 

In order to get the best approximation to the true 
arrival distribution over a 24 hour period on board ship, a 
group of experienced Naval personnel were surveyed for their 
estimates of the actual arrival pattern (i.e. hydraulic 
loading) of water closets and urinals over a 24 hour period. 
From such a survey, an empirical distribution was developed 
by employing a technique similar to that used to get job 
time estimates for PERT programming (Program Evaluation and 
Review Technique). The details of this method of estimation 
and how it was used are contained in Appendix B. The re- 
sults are shown in Figure 4a. 

Figure 4b gives a similar hydraulic loading pattern 
derived for a commercial cargo ship with a 40 man crew 
[Ref. 13, p. 58]. Although the routine on a commercial 
cargo vessel is considerably different from that on a 
military vessel, it can be noted that "peak loadings are 


+ 
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Figure 4a. 
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Figure 4b. 
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Estimation of Distribution of Sewage 
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double the average loading and almost five times the minimum 
loading" for both distributions, and that, in general, they 
compare favorably. 

Once this distribution was established, it was 
necessary to convert it to realistic quantities of sewage 
per hour. Again, the hydraulic shipboard data obtained in 
MEL Report 346/64 [Ref. 19] (which is the basis of current 
Navy design parameters of 30 g.m.p.d.) was used to estimate 
a distribution of the daily per capita sewage flow. This 
data, however, was revised in a thesis by D. B. Campbell 
(a fact not incorporated in design specifications listed in 
Ref. 2, p. E-1). Combining both sets of information it was 
decided that daily per capita sewage flow rates could best 
be represented by a beta distribution with a most likely 
(mode) value of 14.8 gal/man/day, a maximum of 20.0 g.m.p.d. 
and a minimum of 11.3! As shown in Appendix C, these values 
were used to estimate the parameters of a beta distribution 
[Ref. 17, p. 196]. 

Observations of the daily sewage generation rate per 
man were then generated from this beta distribution for each 
day of the simulated week. Once a daily average was gen- 
erated it was converted to an hourly average for a ship 
with 240 men. (The number of men was arbitrary and could 
be read in by the programmer.) Finally, it was translated 
to the empirical distribution by conversion ratios of the 
hourly relative value to the relative mean. For example 


at 1200, if the beta realization was 15 gal/man/day then: 
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Relative Value 
Hourly Realization = — 
Relative Mean 


Daily Beta Variate 





Жошы a а ПО ВОР еа (12) 
24 
d 15 
Realization = х= = 240 = 232 gal 
rM 24 


Such values were the hourly realizations of X(t) used 


as the input to the basic model described in subsection 1. 
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IVA RESULTS 


The results of simulations one and two are presented on 
pages 32 through 35. Figure 5a and 6a are graphs of the 
average maximum level of sewage held as a function of pro- 
cessing rate for the five policies tested by simulations one 
and two reSpectively. Figures 5b and 6b represent the 
average total time the processor must run in a week's period 
(of a possible 168 hour maximum), again as a function of 
processing rate for each of the five policies. The results 
agree well with intuition. 

Looking first at the results common to both simulations, 
one observes that the maximum average Sewage level has an 
inverse relationship to the processing rate. The ОО 
savings in tank size decreases with higher processing rates, 
and, for rates greater than 210 gal/hr, the savings appear 
to be negligible. For policies in which processing time 
can vary, the total weekly processing time is also inversely 
proportional to the processing rate, and, over the range of 
rates considered, the function is nearly linear. 

The different policies considered have a significant 
effect on sewage level and processing time. Policy two 
describes a limit on minimum sewage level obtainable for 
the rates considered. The curve of maximum sewage level 
increases with each higher policy number up to policy five. 
Conversely, for average processing time the trend is just 


the opposite as was expected. Pollcyscwenmrc un p 
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maximum processing time and each succeeding policy number 
requires less. For the policy of strictly holding sewage 
(policy 1), the average sewage accumulated in a week's 
period (assuming the ship to be a destroyer with a comple- 
ment of 240 men) was 21,100 gallons, with a variance of 

8.2 x 10* gallons? for simulation опе, апа 21,080 gallons 
in ce o ita a) ons for simulation two. 
The average values for a 12 hour holding period were 3,057 
gallons for simulation one and 2,700 gallons for simulation 
two. 

Certain additional observations could be made that were 
peculiar only to simulation one. First, a definite change 
in the slope of the curve for maximum tank level vs pro- 
cessing rate occurs at 150 gal/hr, with little additional 
savings in tank size being realized beyond that. For pro- 
cessing rates above 150 gal/hr, the average maximum tank 
level was never above 2,500 gallons. The variance/mean 
ratio for maximum tank level is on the order of from 2 to 10. 
It is usually highest for a processing rate of 140 gal/hr, 
and drops off immediately thereafter. For policy four the 
ratio is much lower than for other policies. For processing 
time vs processing rate the variance/mean ratio is on the 
order of 10°*, and policies 3 and 4 show little appreciable 
difference in time savings. 

For Simulation two, it can be seen that, with policy 
two aS a reference, the curve of maximum level vs processing 


rate levels off more rapidly for policies 3 and 4 beginning 
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at approximately 160 gallon/hr. The curves for policies 3 
and 4 follow the same contour with policy 3 having a maxi- 
mum sewage level on the order of 1000 gallons less than 
policy four for all rates. The variance/mean ratio for 
maximum sewage level is on the order of 10° to 10° at 140 
gal/hr and decreases as processing rate increases. At 
higher..processing.rates.the ratio is significantly less for 
Policy four (approx. 20) than for other policies. For pro- 
cessing time vs processing rate, policies three and four 
have the same slope, with policy four showing only a slight 
Savings in time of about 5 hrs for the week. The variance/ 
mean ratio for either policy is again very small, on the 
erder of 0.5. 

A comparison of the two simulations reveals that, al- 
though the average results of simulation one are on the 
same order as those of simulation two, the variance of the 
results is much smaller for the former. This is not supris- 
ing, since simulation one is based on the Poisson process 
in which the variance of the asymptoptic distribution of 
N(t) will equal the mean. In the case of a non-homogeneous 
process in which the parameters do not differ greatly, it 
can also be expected that as the number of ie in- 
creases the sample variance will approach or value on the 
Same order of magnitude as the mean. For simulation two, 
ehe realizations Of the dally average sewage ssateue ema 


сап varv from 11.3 gal/man/day to 20.0 gal/man/day which is 
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then imposed on an empirical distribution with even greater 
variance. The result is a process with reasonably large 
variance; certainly much greater than that resulting from 


simulation one. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. INTERPRETATION OF RESULTS 

In order to prescribe the most feasible combination of 
processing rates, policies and holding tank size, more 
information regarding the cost function of each particular 
parameter would have to be developed. If costs were as- 
‘=sumed"to follow some -inereasing.linear function.of pro- 
cessing rate, holding tank size and processing time, then 
the policy and rate which would minimize a measure of ef- 
fectiveness taking all such costs into account would be 
expected to be most feasible. 

Based on the results of Chapter IV and with the under- 
standing that cost functions must be applied, it appears 
that a system consisting of a holding tank of 2,700 gallon 
capacity for sewage combined with a processing unit capable 
of handling 160 to 170 gallons of sewage per hour would be 
most beneficial. Further such a system should be operated 
under a processing policy similar to that of policy 3. 

It was contended in Chapter I that the Navy's design 
parameter of 30 gallons/man/day was over estimated. The 
simulations, which were developed using the best data that 
could be obtained, bear out this conclusion. More specifi- 
cally, if this current value were used to construct a 12 
hour holding tank for a destroyer of 240 men, the minimum 
tank size would be 3600 gallons. The values obtained 


through simulation, which were measured during the period 


* 
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when traffic flow was the greatest (0500 to 1700), were, 
again, 3,050 gallons for simulation one, and 2,700 gallons 
for simulation two. Even with the addition of a 200 gallon 
safety factor, it appears evident that a savings of between 
400 to 700 gallons could be realized for each 12 hour 
holding tank ee Even greater savings would be 
„expected..for larger tanks. 

The capacity of a holding tank ren on the quantity of 
diluent (flushing liquid) used can also be computed using 
an equation presented in a report by the Maritime Research 


Information Service [Ref. 13]. Namely: 


С = 1.25 x M x (BAD) | (13) 
where 

C = tank capacity 

M = man days of storage 

B = daily per capita body waste = 0.437 gallons 

D = diluent in gallons/man/day. 


If the value of diluent is taken to be 15 gal/man/day, 
the unmodified mean of the data obtained by MEL Report 
346/64, then the recommended holding tank size would be 
2,320 gallons. The values obtained by simulation would 
thus appear conservative; however, it is important to 


realize that they also account for peak arrival periods. 


B. RECOMMENDATIONS 
The simulations described in this thesis offer sufficient 
merit for use in the design of sewage treatment systems. It 


is recommended that a study Бе паче олово 
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statistics on the arrival pattern of head facilities aboard 

a Naval ship. Such new data could then replace the empirical 
distribution developed for simulation two, or be used to 
revise the intensity function of simulation one. This 

should lead to even more meaningful results than those 
presented. 

It is further recommended that design parameter for 
daily per capita sewage generation be reviewed. If it is 
decided that it is in fact overestimated, then it is recom- 
mended that 12 hour holding tanks of the size presented in 
this report be tested on board destroyer-sized ships, and 
that the simulations be employed to develop holding tanks 


for larger ships as well. 
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APPENDIX A 


CALCULATION OF ARRIVAL PARAMETERS 


As pointed out in the body of the thesis, the mean value 
ШЕШСЕТбп was used to есета се Ебер Бесе аео R x 
1,2,3). This was done by computing average daily and weekly 
sewage generation first by using the average values pre- 
sented in MEL Report 346/64 [Ref. 19], and second by use of 
policy one and simulation two. The values obtained from 
weekly observations of sewage build-up as estimated by 
simulation two were used since 1) they were considered more 
representive and 2) they would also serve as a means of 
normalizing both simulations so that over long periods of 
time the two simulations would give similar results for 
strict holding. Average daily and weekly sewage generation 
was then divided by 4.9 gal/arrival to determine the number 
of arrivals/period. 


The results were: 


Weekly average 21,000 gal/week/4.9 gal/arrival 


4260 arrivals/week 


similarly: 


dally average = 730 arrivals/day. 
These were then equated to the mean value function which 


was 


M(t) = 5(m(t)) + 2 (m(t')) = 5[2, (t,)+42 (t2-t1)+43 (t-t2)] 


> 3: PNP Ges veh) este.) | (А-1) 
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for the weekly value, and 


nit) <a) (A-2) 
for the daily value. | 
Substituting in the values forti where. US 
288 five minute increments and setting them equal to the 


weekly and daily sewage build-ups, one gets: 


weekly: 6604, + 7681, + 588), = 4260 (A-3) 
daily: 132%, + 72%, + 844, = 730 (A-4) 
where: t, = 132, t, = 204, t = 288, t,' = 132. 


In order to solve the above simultaneous equations a 
third equation was needed. This was obtained from weekends 
and computing a third average build-up for that period to 


obtain: 


2m(t') = 2,100/4.9 = 420 arrivals/weekend 


or 


2641, + 3121, = 420. (A-5) 
Solving (A-3) through (A-5) simultaneously one obtains 


the final estimates of À; (1 = 1,2,3) namely: 


^ 

A, = 4.5 arrivals/5 min 
A, = 1.0 arrivals/5 min 
ху = eee el aio акы си 


which agree well with intuition. 
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APPENDIX B 


METHOD OF ESTIMATION FOR AN EMPIRICAL DISTRIBUTION 


In order to obtain the hydraulic loading pattern of 
water closets and urinals aboard Naval ships, it was nec- 
essary to resort to a subjective estimation technique. 

The technique used is the same as that for obtaining job 
time estimates for PERT scheduling computations [Ref. 14]. 
Such a technique is employed when one is unable to do any 
statistical sampling. 

It was, therefore, decided to survey a cross-section 
of experienced Naval personnel, both officers and enlisted, 
for their subjective estimates of the relative usage of 
head facilities for each hour of a 24 hour period. Their 
estimates were made on a relative scale of 0 to 20, (0 in- 
dicating the facility was not used at all, 20 indicating 
it was used to capacity). Each person was to give a high, 
most likely, and low estimate, independent of estimates 
made by others. 

It has been "historically accepted" that job duration 
times described by such estimates are beta distributed over 
a finite interval given by the high and low estimate. The 
same distribtuion can be assumed to apply to estimates of 
the arrival pattern. The expected value of such a distri- 


bution can then be approximated by: 


U, = (at4m+b) /6 | (B-1) 
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and the variance by; 


Wo [ (b-a) /6] * (B-2) 
where: | 
а = low (optimistic) estimate 
m = most likely (mode) estimate 
= high (pessimistic) estimate. 


Using the expected value, U as the best estimate of each 


er 
individual. А11 such values were then averaged to give the 
overall best estimate for each hour of the 24 hour period. 


The results are the empirical distribution plotted in the 


body of the thesis (Figure 4a). 
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АРРЕМОТХ С 


ESTIMATION OF BETA PARAMETERS 


A beta distribution that is transformed to some inten M 


пег һап (0,1) 15 given by: 


Г(А+В) 
——— (b-a)i-A-B (q-a)À-l (p-T)B-1, p>T>a 
Г(А)Г(В) 
f(T) = l (С-1) 
0 elsewhere 


where T is the transformed beta variable given by: 


Т = a+(b-a)X (С=2) 
X being a (0,1) beta variable and a,b being the lower and 
upper bound of the transformed variable. If the high (b), 
low (a) and most likely (m) values of the distribution can 
be determined, then the transformed beta parameters can be 


related to (0,1) beta parameters according to: 


А = 46E(X) [E(X) -E(X)2 - =] (c-3) 
A-AE(X) 
B - ———— (C-4) 
Е(Х) 
where 
4m+b-5a 
E(X) = — (С-5) 
6 (b=a) 


Using the values obtained from MEL Report 346/64 1.e.: 


11.3 gal/man/day 
14.8 gal/man/day 
20.0 gal/man day. 


os 
uU H H 
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Then the values for the parameters can be solved for, using 


the above equations to obtain: 


А = 4.3 В = 5.7 
These values are then used to generate a beta variable, X, 
by first generating two independent gamma variates Y,~G(1,A) 
Баар (l,B) and computing the ratio 

Yi 
X = (С-6) 
YX +X, 

It then remains only to transform X to a realization of T 


using equation (C-2). 
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APPENDIX D 


USE OF SIMULATION ONE AND TWO 


The two simulations were programmed in Fortran IV for 
use on the IBM 360/67 computer. Both simulations allow the 
option of plotting the graph of the sewage build-up (in- 
ventory) for each hour of the week. (Examples of such 
plots are given under computer output.) If one uses the 
option, however, he must insure that the appropriate JCL 
cards are used for his particular hardware. 

Only one data card is needed for either simulation. 
All that must be provided are the parameters to be tested, 
an indication of print routine to be used and a flag for 
calling the plot routine. For each simulation the format 
of the data card is as follows: 

a) Simulation One. The following variables must 


be supplied: 


Variable Cols on Card Description 

Р1ап 1 0-5 As described in the body 
of the thesis, these are 

Plan 2 6-10 the arrival rates/5 min 
that are used for the 

Plan 3 11-15 Poisson generator. 

Alpha 1 16-18 The initial processing 
rate seed, in gal/5 min 
processed. It is in- 


creased by 10 gal/hr as 
soon as it is read in. 


NPLAN 19-20 A two digit integer 


specifying the processing 
plan to be tested. 
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NOPRIN 21-22 A two digit integer to 
indicate which of four 
print routines is to be 
used. (See below for 
discription of these 
routines.) 


NO DRAW 23-24 A flag to indicate if 
plot routine is desired 
(any 2 digit number) or 
not (use zero). 


ZSTOP , ZSTART 25-20 These are the level 
30-34 parameters for policy 4 
їп ЕСКЕ а 
TSTOP 34-36 A 3 digit integer repre- 
senting the stopping 
time for policy 5. 

For simulation two the same values are read in with the 
exception of the Poisson parameters, instead the number of 
men representing the size of ship is read in as a 3 digit 
integer. 

For the print routines it should be noted that if full 
printout is used that this can add considerably to the run- 
ning time of the program. Depending on the read in for 
NOPRIN, the routines are: 

NOPRIN ROUTINE 

1 Print only the final statis: 
tics for ZMAX matrix and the 
processing time matrix. 

2 Print above plus the times 
that the processor stop and 
restarted during the week. 
Print also the processing 
rate, run number and maximum 
tank level on each run. 

3 | Print the above plus the 


entire Z(t) matrix (also the 
x(t) matrix for simulations DE 
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4 Print the above plus the 
X(t) matrix and arrivals 
matrix. (Simulation One) 

In conclusion, it should be pointed out that simulation 
one requires approximately 2000 random numbers to be gen- 
erated per replication; it requires considerably more com- 
puter time than simulation two. This is inpart compensated 


for, however, by its smaller variance, so that fewer runs 


are needed. 
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SUBROUTINE RANDU(IX, IY,YFL) 


un GENERATES RANDOM NUMBERS FOR USE IN THE MAIN 
COMMON I,J,ITIME,NOPRIN 
IY= I X#&5539 
IF(1Y)5,6,6 
IY=1Y+2147483647+1 
YEL=IY 
YFL=YFL*.4656613E-9 
RETURN 
END 
Comment cards denoted by and including * should be 


removed prior to running programs. 
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